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چکیده 


تخمین بار معلق در طیف وسیعی از مسائل از قبیل طراحی مخازن سدهاء برآورد میزان فرسایش و رسوب‌گذاری اطراف پایه‌های پل و مدیریت 
حوضه‌های آبخیز به‌کار گرفته می‌شود. در این پژوهش به‌منظور تخمین بار معلق حوضه دره‌رود. مقادیر دبی و بار معلق در ۱۶ ایستگاه هیدرومتری طی 
دوره مشترک آماری ۱۵ ساله (۱۳۸۰-۹۴) جمع‌آوری گردید. پنج الگوی مختلف بر اساس میزان تأثیرگذاری متفیرهای دبی و خصوصیات فیزیوگرافی 
زیرحوضه‌ها شامل مساحت شیب ضریب شکل و شماره منحنی بر بار معلق حوضه تعریف شد. ضمناً با در نظر گرفتن پارامترهای مساحت و شیب 
زیرحوضه‌ها به دو گروه اول و دوم تقسیم‌بندی شدند. عملکرد مدل‌های شبکه عصبی مصنوعی (۸» سامانه استنتاجی فازی- عصبی تطبیقی 
(۸(3۲15) و برنامه‌ریزی بیان ژن (0101)) در پیش‌بینی بار معلق مورد بررسی قرار گرفت. نتایج نشان داد تخمین بار معلق با به کارگیری الگوی ترکیبی 
شامل کلیه خصوصیات فیزیوگرافی و دبی با پیشترین دقت همراه بود. در بین مدل‌های هوشمند بهترین عملکرد متعلق به مدل "00101 بود. در گروه اول 
این مدل بیشترین ضریب تعیین (۰/۶۸-**41 کمترین مقدار ريشه میانگین مربعات خطا ((0/02) 0۷15]7۷/۶۹) و بیشترین ضریب نش-ساتکلیف 
(۹-۰/۵۵) را در مقایسه با سایر مدل‌ها به خود اختصاص داد. در خصوص گروه دوم نیز مدل 0۳۳ با دارا بودن مقادیر 182 ۷5۳( و 5 به‌ترتیب 
برایر با ۰۰/۷۲ ۹۷۵/۲۶ و ۰/۴۳ برتری محسوسی داشت. با استفاده از مدل 0 برای گروه‌های اول و دوم مدل‌های منطقه‌ای رسوب استخراج شد. 
طیق کا ط ال ها ۱۲۸۸۰۵۱۴ الک راو منکیمن ۶۱۳۲ لبون گس رات سای توس شبکه آیزاهه ها خر کل بحرشه دزم روف خاوتها شیور 


به‌طور متوسط سهم هر کیلومتر مربع حوضه, حدود ۱۰۰۰ تن بوده است. 


واژه‌های کلیدی: بار معلق. حوضه دره‌روده خصوصیات فیزیوگرافی» دبی ویژه رسوب. مدل‌های هوشمند 


مقدمه 


در رودخانه‌ها همواره فرسایش و انتقال رسوب انجام می‌گیرد. از 
این‌رو بررسی ظرفیت حمل رسوب جریان و سازوکار انتقال رسوب در 
هیدرولیک رودخانه‌ها و ریخت‌شناسی آن از اهمیت قابل توجهی 
برخوردار است. انتقال درات رسوبی (بار کل) در مجاری روباز به دو 
شکل بار بستر و بار معلق صورت می‌گیرد. بار بستر بخشی از بار کل 
خروجی حوضه آبخیز بوده که در تماس با بستر حرکت می‌کند و بار 
معلق بخشی از بار کل بوده که به‌علت آشفتگی جریان به‌صورت 


دانشکده کشاورزی و منابع طبیعیء دانشگاه محقق اردبیلی 

(#- نویسنده مسئول: ۱۵۴ ۵ ۲نا 6۵0ص _قط :لتمصط) 

۴- دانشیار دانشکده کشاورزی و منابع طبیعی مغان. دانشگاه محقق اردبیلی 
1 1022067/5۷ :]10 


غوطه‌ور همراه با جریان آب حرکت می‌کند (۱۵). امروزه پدیده انتقال 
رسوب به‌عنوان یکی از بزرگ‌ترین مشکلات بهره‌برداری از منابع آب 
های سطحی مطرح می‌گردد. به‌طوری‌که باعث کاهش ظرفیت 
آبگذری آبراهه‌ها و تاسیسات انتقال آب و هم‌چنین تغییر کیفیت آب 
به لحاظ مصارف شرب و کشاورزی می‌شود (ع)» به دلیل تأثیر 
پارامترهای مختلف بر انتقال رسوبات در رودخانه‌هاء روابط تجربی 
متعددی در خصوص پیش‌بینی بار معلق توسط متخصصان ارائه شده 
اما در بسیاری از موارده نتایج حاصل از روابط تجربی با مقادیر واقعی 
مطابقت نداشته و لذا هیچ کدام مورد مقبولیت جهانی قرار نگرفته است 
(۵). امروزه استفاده از مدل‌های هوش مصنوعی به‌عنوان راه‌کاری 
جدید برای تحلیل و درک رفتار غیر خطی مسائل متنوع مهندسی. 
ازجمله پدیده‌های هیدرولیکی و هیدرولوژیکی گسترش يافته است. در 
زمینه تخمین بار معلق با استفاده از مدل‌های هوشمند پژوهش‌های 
فراوانی صورت گرفته است که به چند مورد از آن‌ها طی سال‌های 
اخیر اشاره می‌شود. شیخعلی‌پور و همکاران (۲۴) از مدل‌های شبکه 


۸ نشریه آب و خاک. جلد ۰۲۴ شماره ۰۴ مهر - آبان ۱۳۹۹ 


عصبی مصنوعی! (۸۲۷(۷)» سامانه اتنتاجی فازی- عصبی تطبیقی " 
(۸(۷۳15) و برنامه‌ریزی بیان ژن " (015۳) جهت برآورد با معلق 
رودخانه سیستان استفاده کردند. ارزیابی دقت مدل‌ها با توجه به 
شاخص‌های آماری خطا نشان داد که هر سه روش نتایج مناسبی در 
برآورد بار رسوبی ارائه می‌دهد و مقایسه عملکرد مدل‌ها حاکی از 
برتری نسبی مدل 13۳) نسبت به سایر مدل‌ها بود. قربانی و دهقانی 
(۱۱) کارایی مدل‌های ((۸ و شبکه عصبی بیزین" (((3) را برای 
تخمین رسوبات معلق رودخانه‌ها مورد ارزیابی قرار دادند. ایشان از 
داده‌های رودخانه سیمینه‌رود واقع در آذربایجان غربی طی دوره آماری 
(۱۳۵۴-۸۳) استفاده کردند. ننایج تحقیق آن‌ها حاکی از توانایی بالای 
مدل 13 در تخمین بار معلق و هم‌چنین مقادیر کمینه و بیشینه 
رسوب بود. نیک‌پور و ثانی‌خانی )۱٩(‏ برای پیش‌بینی رسوبات ایستگاه 
هیدرومتری مشیران بر روی رودخانه دره‌رود از مدل‌های 0۲۲ 
5 و ٩۷‏ استفاده کردند. نتایج پژوهش آن‌ها بیانگر عملکرد 
قابل قبول مدل‌ها و برتری مدل 5۷8 در مرحله صحت‌سنجی بود. 
ابوالفتحی و همکاران (۱) میزان رسوب رودخانه گلرود واقع در استان 
لرستان را با استفاده از انواع مدل‌های ((۸ شامل تابع پایه شعاعی* 


381 و پرسپترون چند لایه " (0۸1) و هم‌چنین رگرسیون خطی 
چند متغیره مورد بررسی قرار دادند. آن‌ها علاوه بر مقادیر دبی جریان 
از ویژگی‌های مورفولوژیکی حوضه شامل مساحت» طول, محیط و 
شیب رودخانه نیز برای مدل‌سازی استفاده کردند. در تحقیق مذکور 
مدل ]۷۲ از تطابق بهتری با مقادیر مشاهداتی برخوردار بود. اسدی 


و فتح‌زاده (۳) توسط شش مدل هوشمند شامل "1 نزدیک‌ترین 
همسایه» مدل درختی ۷5 فرآیند گوسی" (0۳)» ماشین‌بردار 
پشتیبان* (45۷14 ماشین بردار پشتیبان تکاملی " (25۷۷) و شبکه 
عصبی پس‌انتشار خطا (۱3۳) مقادیر بار معطق هشت ایستگاه 


هیدرومتری واقع در استان گیلان را برآورد نمودند. ارزیابی دقت مدل- 
ها با توجه به شاخص‌های آماری خطا نشان داد که مدل فرآیند 


گوسی در مقایسه با سار مدل‌ها از کارآیی بهتری برخوردار بود. دولت 
کردستانی و همکاران (۷) برای پیش‌بینی بار معلق روزانه ایستگاه گرو 


1- ۸۱1110121 ۲6۵1 ۱۱60۱۷0۲1 )۵۵۵( 

(۸۵۲۳]۵) اور 5 18/۵۲۲۵۲۵۵6 ۲۵۲۵-۲۵22 ۸020۷۵ -2 
(طنان) عمتصصهتعن۲۳۵ مملووع۲)۵۲ ععع) -3 

(۲ظ) ۱۱6۱0۲1 لویه۲ مهاوعظ -4 

(۳ظ۳) ممتاعصی۲ متمظ ۴۵0121 -5 

(صا۱) مماو۲06عظ ررض آ-تالن۱۷۲ -6 

۲صطونه۱ )۱6۵769 ک -7 

8- )120959128 5 

(5۱۷۱) مصتط۷]۵ ۷۵۵۲۵۲ 5۱۵۵0۲۲ -9 

فصتطمقظ ۷6۵۵۲ ۹0۵۵۵۵۲۲ ۷۵11210۳۲۵۲۲ -10 


واقع در استان هرمزگان از مدل‌های ۸ 9 استفاده 
کردند. در تحقیق آن‌ها با در نظر گرفتن تأخیرهای دبی و رسوب سه 
الگوی ورودی مختلف تعریف گردید. نتایج تحقیق آن‌ها نشان داد که 
مدل ۸۱۴15 به‌ازای هر سه الگو از عملکرد بهتری در برآورد بار 
مدل بهینه برای برآورد میزان بارمعلق رودخانه هلیل‌رود براساس دبی 
جریان به مقایسه مدل‌های منحنی سنجه رسوب. ۸۲۲ و ۸۳۲15 
پرداختند. نتیجه کلی حاصل از تحقیق مذکور مدل ۸۳۲ را به 
عنوان بهترین روش برای شبیه‌سازی میزان رسوب معلق معرفی نمود. 
صمدزاده و همکاران (۲۱/ برای ایستگاه‌های هیدرومتری واقع در 
(۱۶) کارایی سه مدل هوشمند شامل ۵۲۲ ۸۳۱۳۲5 و ۳۲ را در 
بررسی قرار دادند. آن‌ها برای مدل‌سازی از مقادیر بارش و دبی روزانه 
نیز استفاده کردند. نتیحه تحقیق آن‌ها حاکی از عملکرد بهتر مدل 
در مقایسه با دوم دل دیگر بود. ایلماز و همکاران (۲۶) برای 
برآورد بارمعلق رودخانه‌ای واقع در ترکیه از مدل‌های "۷1۸1۲5۲ و 
کلونی زنبور عسل ۲ استفاده کردند. آن‌ها با به‌کارگیری داده‌های دو 
برای آموزش مدل‌ها تعریف کردند و در نهایت مدل ۱1۸185 به‌عنوان 
دقیق‌ترین روش شناخته شد. سأمانتارای و قوس (۲۳) از سه نوع مدل 
معلق حوضه آبخیزی واقع در هندوستان استفاده کردند. مقایسه نتایج 
نشان داد که مدل ۱۷۱۲ با دارا بودن بیش‌ترین ضریب 1 و کم‌ترین 
۳ از بهترین عملکرد برخوردار بود. حسن‌پور و همک‌اران (۱۲) 
برای تخمین بار معلق رودخانه سیستان از مدل‌های منحنی سنحه 
رسوب. ۵ ٩۷۱۷‏ و 70۷-۹۷۷۲۲ استفاده کردند. با در نظر 
گرفتن شاخص‌های آماری مختلف» برتری مدل هیبرید ۳۷۲-۹۷۷۲ 
نسبت به سایر مدل‌ها نمایان بود. صالح و همکاران (۲۰) با به 
کارگیری مدل‌های ۷510116 ۷۲5۳ و 27 اقدام به شبیه‌سازی 
براساس مقادیر دبی جریان و بار معلق الگوهای مختلف تعریف گردید. 
نتایج تحقیق آن‌ها نشان داد که هر سه مدل از دقت قابل قبول برای 


(رفخها) مصتاموگ مملووءتعومف۴ ۸02۵۷۵ ماهلیه ت۷۲ -11 
7 866 -12 
(۳۵۷۲-۹۷) ۷۱۷۲ 5-مصتماوناه صقعصصت) ۲۱22۷ -13 


جزئی بود. عدنان و همکاران (۲) برای پیش‌بینی بار معلق رودخانه 
جیالینگ واقع در کشور چین از مدل‌های ۲۳۵۲۳۲ و -۸(۳۲5 
۷ استفاده کردند. برای این منظور از داده‌های دو ایستگاه 
هیدرومتری در حوضه آبریز مذ کور استفاده شد. نتایج حاکی از آن بود 
که مدل 1۳۳۲5 می‌تواند با دقت بالایی برای براورد بار معلق هر 
دو ایستگاه به‌ کار گرفته شود. 

حوضه آبخیز دره‌رود یکی از مهم‌ترین حوضه‌های استان اردبیبل 
بوده و رودخانه‌های این حوضه تأمین کننده آب شرب اغلب شهرهای 
استان (اردبیل نمین» نیره سرعین و مشگین‌شهر) می‌باشد. علاوه بر 
اين. دو سد مخزنی بزرگ (سبلان و یامچی) در حال بهره‌برداری و 
چند سد مخزنی در حال ساخت در سطح حوضه وجود دارد. نظر به 
موقعیت استراتژ 
برداری از آن. مدل‌سازی و پیش‌بینی میزان رسوبات معلق در این 
حوضه می‌تواند در راستای مدیریت منابع آب استان اردبیل نقش به 
هوشمند شامل ۸۸ ۸۳۲5 و 0 در تخمین بار معطق کل 
حوضه دره‌رود با به‌کارگیری خصوصیات فیزیوگرافی و دبی رودخانه‌هاء 
مورد ارزیابی قرار گرفت. 


مواد و روش‌ها 
منطقه مورد مطالعه 

استان اردبیل دارای حوضه‌های آبخیز دره‌رود بال‌هاروده قزل‌اوزن 
و ارس می‌باشد. حوضه آبخیز دره‌رود با مساحتی بالغ بر ۱۴۳۱۰ 
کیلومتر مربع» بزرگ‌ترین حوضه استان بوده که ۷۶ رودخانه دایمی و 
فصلی به طول ۲۰۶۷ کیلومتر در آن جریان دارد. ۶۸ درصد وسعت 
شرقی است. مهم‌ترین رودخانه‌های این حوضه عبارتند از: ارس» دره 
رود بالخلی‌چای» قوری‌چای و قره‌سو. در پژوهش حاضر به‌منظور 
تخمین بار معلق کل حوضه دره‌رود از مقادیر اندازه‌گیری شده دبی و 
بار معلق در ۱۶ ایستگاه هیدرومتری واقع در سطح حوضه مذکور 
استفاده شد. داده‌های دبی و بار معلق از شرکت آب منطقه‌ای استان 
اردبیل اخذ شد. شکل ۱ موقعیت جفرافیایی ایستگاه‌های هیدرومتری 
واقع در حوضه دره‌رود و شکل ۲ نقشه آبراهه‌های این حوضه را نشان 
می‌دهد. 
داده‌های مورد استفاده 

اندازه‌گیری بار معلق در ایستگاه‌های هیدرومتری منطقه مورد 
مقادیر دبی جریان به‌صورت روزانه ثبت شده بود. در این پژوهش» 
دوره آمارشن ثبت شده طی ۱۵ سال اخیر (۱۳۸۰-۹۴) به‌عنوان دوره 


کاربرد سه مدل هوشمند در بر آورد بار معلق حوضه‌های آبخیز ... ۸۲۹ 


مشترک آماری انتخاب گردید. بنابراین پس از حذف داده‌هایی که 
خارج از بازه زمانی مشترک بودند» مجموع تعداد مقادیر مشاهداتی بار 
معلق به‌ازای کل ایستگاه‌ها برابر با ۲۱۲۱ داده بود. لازم به ذکر است 
که اندازه‌گیری بار معلق توسط نمونه‌بردار دستی 1۳1-48 انجام گرفته 
بود (شکل ۳). در جدول ۱ خصوصیات فیزیوگرافی زیرحوضه‌های 
بالادست ایستگاه‌های هیدرومتری حوضه دره‌رود شامل مساحت (ح» 
ضریب شکل (,5) ارتفاع متوسط از سطح دریا (4۳ زمان تمرکز (:4 
شیب آبراهه اصلی (8) و شماره منحنی (()) گزارش شده است. 
برای تهیه نقشه کاربری اراضی ابتدا با توجه به دانش موجود درباره 
پوشش اراضی منطقه و روش تفسیر چشمی تصاویر 0ات2 00216 
هفت طبقه کاربری اراضی شامل مرتع خوب مرتع متوسط و ضعیف» 
و باغ در منطقه مورد مطالعه تعیین گردید. در مرحله بعده در محیط 
نرم‌افزار ۸6615 نقشه گروه هیدرولوژیکی و نقشه کاربری اراضی با 
هم تلفیق شد. سپس با استفاده از رابطه ۱ نقشه شماره منحنی حوضه 
با توجه به جدول استاندارد ۹6٩‏ به‌دست آمد. 

)۱( (6۸ <) < < ۲7 
درصد مساحتی از حوضه است که شماره منحنی آن :) باشد. برای 
محاسبه زمان تمرکز حوضه در محیط نرم‌افزار ۷۷5 ابتدا زمان 
تأخیر توسط رابطه سازمان حفاظت خاک آمریکا (505) محاسبه شد 
(رابطه ۲) 


1(7+) نا 
(۲ از 
1900/3 4 


در رابطه ۲ یر" زمان تأخیر حوضه (ساعت) .1 طول رودخانه 
اصلی (فوت) ا شیب متوسط حوضه (/) و 5 نمایه نگه‌داشت رطوبت 
توسط خاک در داخل حوضه (اینچ) است که مقدار آن برابر است با: 
۳ 4 - سظ < و 

پس از محاسبه زمان تأخیر زمان تمرکز حوضه (:1) توسط 
معادله (۴) به‌دست آمد: 
9 وم *< 1.67 < .1 

به‌منظور تسهیل مدل‌سازی بار معلق حوضه دره‌رود و رسیدن به 
نتایج بهتر ایستگاه‌های هیدرومتری این حوضه به دو گروه تقسیم 
بندی شدند. برای تقسیم‌بندی ایستگاه‌ها دو پارامتر مساحت تحت 


پوشش و شیب زیرحوضه‌های بالادست ایستگاه‌های هیدرومتری مد 
نظر قرار گرفت. به‌طوریکه زیرحوضه‌های کوچک با شیب زیاد 
(بیشتر از ۱) که متعلق به سرشاخه‌های حوضه دره‌رود است در گروه 
اول و زیرحوضه‌های متوسط و بزرگ با شیب کم (کوچک‌تر مساوی 
۱ در گروه دوم قرار گرفتند. طبقه‌بندی مربوط به گروه‌های اول و 
دوم در جدول ۱ ارائه شده است. 
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47 42 4.2 48.048768 48.358336 -.-4 
۷ 


2-308 


239-32908 


عباد۱۷ 


2 ۲:۵ : 6 


۰ 1۵۷ : 8 


ارتفاع 
معا 


"2-2 2-54 
"2-2-23-13594 


دوست بیگلو 
اکتا[ 
۱212100 


2-۹28 
ِ 
3 
156: 
"29 28 


راهنما 
0 1( 


2046 


28 0296 


ایستگاه هیدرومتری + 
«متاهای عتعا۲۵۳۵ ۲3 
مرز حوضه دره رود | __] 
فوم(- ومد ۵۶ معقمتم8 


هم 
0 30 20 510 0 قِ 
۶ ۱625 سا یت آت 
6 
سس 
‌‌ 


- 72 1-392 48048768 ۰ 48398336 ۰ / 4 


شکل ۱- موقعیت جغرافیایی ایستگاه‌های هیدرومتری در حوضه آبخیز دره‌رود 
۳۷2۵0 ۱2۲۵۵-10۱0۷ عمط رز فصمتاهای امصصم بل بو ۵ حمتایم1۱ افنطام‌منو۵ع -1 ۲12۱۳۲۵ 


47 2 4.2 48.048768 48.358336  -.-4 


۷" 
1 89 
۸ 1 
نع آح 
بد 
۹ 9 
[ 8 
ل 
4 
9 9 
[ ِ 
زره 9 
۳2 9 
9 9 
8 5 
8 تست 8 
3 مها ات 
5 نته ین گر ز 1ج 
وو| اعموطه‌هل۸ > 5 
آلادیزگه ۲ 
( 
ایب 0/۱ 
3 
ها ۱۱۳۳۹5 1 
5 / 8 
بدا كِ "1 ِ 
۳۴ 3 
ه‌ / 7 ‌ 
7 مج وم ی ممتاهاه عتتاعصهع1۳۲ 
۲ 
و ۲ رودخانه 
مرز حوضه دره رود ۳۹ 
1 ئ‌ 
30 20 و5 ۵ مه نع ۵۶ توعقمتاهظ ف 
0۱06۲5 سا سس ۵اه ات 
‌ِ 
‌ 
0 


4732 1392 48048768 48:358336 ۰ ۳-4 


شکل ۲- نقشه آبراهه‌های حوضه آبخیز دره‌رود 
0 ۱۵۵-10۱۱ عط) صز منته۲ ۱ ۵۶ ۱۷۲۵ -2 ۲12۱۲۵ 


کاربرد سه مدل هوشمند در بر آورد بار معلق حوضه‌های آبخیز ... ۸۳۱ 


شکل ۳- نمونه‌بردار دستی ۳1-48( برای اندازه‌گیری بار معلق 
0 هداز م6 متس همه ۲۵۲ ۱۱-48 مصتام‌صقی امیصها۱۷ -3 مساعز۲ 


تعداد داده‌های دبی و بار معلق مورد استفاده برای گروه‌های اول و 
دوم به‌ترتیب برابر با ۸۷۰ و ۱۲۵۱ بود. در پژوهش حاضر ۸۰ درصد 
داده‌ها برای آموزش" و ۲۰ درصد برای آزمون" مدل‌های هوشمند به 
کار گرفته شد. لازم به ذکر است که انتخاب داده‌ها برای مراحل 
آموزش و آزمون به‌صورت تصادفی انجام گرفت. در جدول ۲ 
پارامترهای آماری مربوط به داده‌های آموزش و آزمون مقادیر دبی و 
بار معلقق در گروه‌های اول و دوم گزارش شده است. 


شبیکه‌های عصبی مصنوعی (۸۸) 

شبکه‌های عصبی مصنوعی (((۸) یکی از مدل‌های محاسباتی 
نرم " بوده که با هدف یادگیری و دستیابی به دانش نهفته در داده‌هاء 
به‌طور گسترده توسط متخصصان حوزه‌های مختلف علوم مهندسی 
به‌کار گرفته می‌شود. این شبکه‌ها قادر هستند با انتخاب مناسب لایه 
ها و پردازشگرهایی به نام نرون» نگاشت غیرخطی بین ورودی‌ها و 
خروجی‌ها ارائه دهند (۲۵). ساختار معمول مدل ۸۷ شامل سه لایه 
ورودی» پنهان و خروجی می‌باشد. لایه ورودی وظیفه توزیع داده‌ها در 
شبکه را به عهده دارد. در لایه پنهان» داده‌ها پردازش شده و در حین 
آموزش شبکه با استفاده از اطلاعات حاصل از داده‌هاء وزن‌های شبکه 
تعیین می‌شوند. لایه خروجی نیز نتایج را به‌ازای ترکیب ورودی‌ها 
استخراج می‌کند (۲۲). در پژوهش حاضر برای تشکیل مدل ۸ 
ساختار پرسپترون چند لایه (۱11 با یک لایه پنهان در نظر گرفته 


210 -1 
2-1 
عصتانام‌ججمی ]501 -3 


شد. هم‌چنین برای آموزش شبکه, الگوریتم لونبرگ- مارکواردت" به 
کار گرفته شد. به‌منظور آموزش بهتر و سریع‌تر مدل, داده‌ها نرمال 
سازی شدند. در پژوهش حاضر برای نرمال‌سازی داده‌ها از رابطه ۵ 
استفاده شد (۲۷). در اين رابطه ,26 مقدار نرمال شده ورودی :26 
مس و «نسک به ترتیب مقادیر بیشینه و کمینه داده‌ها می‌باشد. 

(۵) ی ال +014 2 ,1 
سامانه استنتاجی فازی- عصبی تطبیقی (۸(۳۲5) 

ژانگ (۱۳) برای اولین بار با به کارگیری ویژگی‌های منطق فازی 
و مدل ۸۸ سامانه استنتاجی فازی- عصبی تطبیقی (۵۳۲5) را 
ارائه داد. در این سامانه, از منطق فازی با هدف بهبود آموزش شبکه 
های عصبی و افزایش قدرت تصمیم‌گیری در شرایط عدم قطعیت 
استفاده می‌شود. مطابق شکل ۴ شش لایه مختلف برای مدل 
5 قابل شناسایی است (۱۳): 

لایه اول (ورودی): شامل متغیرهای ورودی ۸1 و ۸62 می‌باشد. 


(۶) 21 << مررمشام0 

2 < مررمتام0 

لایه دوم (فازی سازی): در این لایه. درجه عضویت () برای هر 
یک از ورودی‌ها بر اساس نوع تابع عضویت در نظر گرفته می‌شود. 

۳ 1,2 < ا,(21) یره < زرریکاجه0 

1,2 ت ,(2 )زوا <ِ زرر‌لامش0 


4- 1,۵۷ ۵۵6۲2- 1 
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شماره 


9 
۱۱۱۵ 


65261 


329 
(۳-114 


1306 
10342 


68928 


6344 
۳.192 


۳-8 


(300 


(۳-144 


13909 


1-19 


339 


رد۱۳۳ 


1-6 


جدول ۱- خصوصیات فیزیو گرافی زیر حوضه‌های بالادست ایستگاه‌های هیدرومتری حوضه دردرود 
۵۹۵ ۱۵۱-1۸0۱۱ عط) ما مصمتاهای عتاء‌صصمنیل بوظ 0 مصتقعوها-طایی رورت م‌تاوزها فص »نطامرهعمتوودام 16 -1 12016 


زمان تمرکز 
ار( 
(بط) متا 


1388 


0234 
68 


139 


60054 


209 


2917 


13.8 


245 


32.972 


1۱-330 


2-2-14 


24 


2-935 


2.16 


2-7 


ارتفاع 
ها( 
)۳ 
14102 


2116 
1112 


1898 


1338 


2227 


1935 


234 


2179 


2117 


1149 


1157 


135 


16085 


166 


1198 


ضریب 


6 
133010 


0۳206 


0295 
033 


(369 


0.579 


۳۹29 


0۳237 


0۳2617 


00201 


038 


0۳46 


0997 


0۳25 


057 


)۳ #6 


(«34 


ِ‌ 


910۳06۰ 
(0 


117 


1.212 
237 


314 
2.98 


2041 


1 
6342 


33 


113 
1031 
(«994 


0۳ 217 


0۳۹13 


0393 


0385 


لایه سوم (قوانین): در گره‌های این لایه» سیگنال‌های ورودی در 
هم ضرب شده و حاصل آن که بینگر قوه انتعال" هر قانون فازی 
می‌باشد. به لایه بعد انتقال داده می‌شود. در رابطه ۴ :۷۷ قوه 


۳ 


مساحت 


۰۹9 
(صعا) 


1.33 


2975 
28 


322-112 


2319 


131202 


نز۱۳۳ 


1۹10029 


16933933 


-<2-62386 


0۳۹402 


06-3929 


4742. 4 


ال۳/ 


2200024 


1--22(9344 


اشتعال قانون .2 ام می‌باشد. 
۳ ِ آ ,(۲2) زولر ۴ (1 نملا < ۷ 


۸) 


نام 
ایستگاه 
ماما 
10 
ننه کر آن 
۱۳۹۵ 
۱۹ 
آتشگاه 
طحعطدد‌اض 
آلادیز گه 
جع 12012 
عنبران 
1( 


ایریل 
111 


هیر 
۳ 
نوران 
۱9۱ 
پل سلطانی 
تصعتاه5 
26 
۱1۹ 
یامچی 
تمصع ۷ 
کوزه تپراقی 
> 
001۳22 
سامیان 
صحتصه 
ارباب کندی 
0( 
ت0صع1 
دوست بیگلو 
اودا0([ 
01200 
مشیرا آن 
صحتتطمه]۷ 
بران 
۱ 


نام رودخانه ‏ رایف 
۱ 10 
6 101۷۵۲ 
۰( 
2۵5-27 ۱۱2۲( 
نوران‌چا ی 2 
ساملا 
قره‌سو 3 
:0 
عنبران 4 
مجتهصا ۸ 
سقزچی‌چای 
-200621 3 
9:7 
۳۱۱۱-۲ 
نوران‌چا ی 7 
/ ۱۹۹۹ 
خیاوچا ی 8 
2۲- 1۳12۷ 
نیرچای 9 
۱۱۱-2( 
بالخلی چا 
بالخلی چای 10 
207 2110111ظ 
قوری‌چای 1 
00 
قره‌سو 2 
:0 
قره‌سو 3 
:0 
قره‌سو 4 
02120-0 
دره‌رود 5 
21۵۳0-۳00[ 
دره‌رود 6 
12700-0] 


نام گروه 
متام 
12106 


گروه اول 
ات 16 
2100۴ 


گروه دوم 
196 
20000 


2100۴۱ 


لایه چهارم (نرمال‌سازی): اين لایه پس از دریافت هر یک از 
قوانین لایه سوم. مقادیر نرمال شده قوانین را توسط رابطه ٩‏ محاسبه 


ی خن 


کاربرد سه مدل هوشمند در بر آورد بار معلق حوضه‌های آبخیز ... ۸۳۳ 


جدول ۲- پارامترهای آماری داده‌های آموزش و آزمون دبی و بار معلق برای گروه‌های اول و دوم 
ص20 6600و 20 )عیق مطه ده 0جه۱ 060صومرویای صرح معبفمل ون 10۲ 0262 اجع) 290 صتقه 0۲ وما200۵ وس آی1ادناه)ه -2 12016 


چولکی انحراف معیار کمینه بیشینه 
دورد لا 50 صصتص]۷ ص۷۱2 
19 1333 0001 1۳10 
2.059 199225 0001 9550.03 
2-332 685 )0 0001 4.450 
6939 1333 0001 123.970 
2259 6369 0001 633100 
113162 ۰85 - 0003 97.912(,-2 
2225 0196 0327 329.00 
2۰394 1-11 0306 2-۳397" 
2 21 


)6 ری 


21 2 


(۵8۵ 


تفسیی‌بندی 

میانگین نوع داده داده‌ها 

۱۱/۳ 1۵۱۵ 
0۵ 


نام گروه 
فا 


دبی 
0500 ۱ 
() داده‌های آموزش 
بار معلق 2 ۲210188 ]1 
15.44 0( 


(00/02) 1020 گروه اول 


ات 16 


تبی 2100۴0 


0305 ۱۱ 
(و) داده‌های آزمون 
بار معلق 12 650[ 
3.350 90( 
(00/020) 1020 
دبی 
8 ۱ 
(و/۳2) 
ار معلق 
0( 
1020 گروه دوم 
(10۳0/020) 0 6[ 
دبی 21000 
5-62 ۱ 
(و) داده‌های آزمون 


داده‌های آموزش 
ها عطصطنه] 1 
999 


بار معلق 12 650[ 


(90 
1020 )00/02( 


"-<2-0917 


شکل - ساختار مدل ۸۳15 (۱۳) 
(13) ۱۵061 ۸۸۱۱۲۲۵ 6( ۵ مساو -4 مساعز۲ 


طبق این رابطه نسبت قوه اشتعال [ امین قانون به مجموع قوه 
های اشتعال همه قوانین محاسبه می‌شود. 
)0 12 < ا, ِ ِِ ۳۷ ِ تروم‌الاماه0 


۷ 


لایه پنجم (غیرفازی سازی ): در اين لایه میزان مشارکت (وزن) 


۱۳ 0 


امین قانون فازی: با حاصل‌ضرب قوه اشتعال" ترمال شده (خروجی 
لایه چهارم) و تابع برآیند (زتبروب+زح) محاسبه می‌شود. 
(۱۰ (۲ + یو + رم۷۷ ع :۷۷ < رتم0 


طاعصمه عط:۲1 -2 
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لابه ششم (جمع‌بندی): در این لابه با جمع مقادیر غیرفازی شدهء 
خروجی مدل ۸۷115 تولید می‌شود. 
(۱۱) ( ۲« :22۷ ع زریاام0 


روش‌های تقسیم‌بندی (افراز) شبکه" و خوشه‌بندی تفریقی " دو 
روش معمول برای طبقه‌بندی داده‌ها در مدل ۸۸۲۳۲5 می‌باشد. افراز 
شبکه فن دسته‌بندی است که در آن هر داده با درجه‌ای مشخص» به 
یک دسته تعلق می‌گیرد. این درجه با درجه عضویت مشخص می‌شود 
(۲۲). تعداد قوانین فازی با افزايش متغیرهای ورودی به صورت توانی 
افزايش می‌یابد و به‌نظر می‌رسد که کاربرد روش افراز شبکه به آزای 
تعداد بسیار زیاد قوانین فازی با مشکل مواجه می‌شود (۱۸). برای به 
کارگیری این مدل توصیه شده است که تعداد متغیرهای ورودی کمتر 
از ۶ باشد (۱۴). با توجه به اين که در تحقیق حاضر تعداد متفیرهای 
ورودی» حداکثر برابر ۵ انتخاب شده. بنابراین استفاده از ایين روش 
جهت پیش‌بینی میزان بار معلق قابل قبول است. 


برنامه‌ریزی بیان ژن (21)) 

مدل 012۳ ترکیب و تعمیم يافته الگوریتم ژنتیک (۵۸) و 
برنامه‌ریزی ژنتیک (0۳) است که بر مبنای نظریه تکاملی داروین " بنا 
نهاده شده است و اولین بار توسط فریرا در سال ۱۹۹۹ ارائه شد (۸). 
در این مدل» کروموزوم‌های خطی با طول ابت و قابل بیان به شکل 
ساختارهای درختی هستند. ساختار خطی کروموزوم‌ها موجب می‌شود 
استفاده از عملگرهای ژننیکی نظیر تکثیره جهش و ترکیب همیشه 
ساختارهای درختی صحیح و معتبری تولید کند به‌طوری که تنها ژنوم 
اصلاح شده به نسل بعد منتقل می‌شود .)٩(‏ لازم به ذکر است که 
ساختارهای درختی از توابع ریاضی و مثلشاتی» عملگرهای ریاضی و 
ترمینال‌ها (متفیرهای مسأله و اعداد ثابت) تشکیل می‌شود. طی فرآیند 
شبیه‌سازی در مدل 6117۲ مراحل ذیل دنبال می‌شود: گام اول» تولید 
جمعیت اولیه از راه حل‌ها است که به کمک متغیرهای مسأله» توابع و 
ترمینال‌ها صورت می‌گیرد (تولید تصادفی کروموزوم‌ها). در گام دوم. 
کروموزوم‌ها به‌صورت بیان درختی ارائه می‌شوند (انتخاب معماری 
کروموزوم‌ها). در گام سوم. کارایی هر عضو جمعیت کروم وزوم‌ها به 
کمک تابع برازش " ارزیابی می‌شود. علت به‌کارگیری تابع برازش, 
حذف جواب‌های ضعیف‌تر و حفظ جواب‌های برتر (نخبه‌گزینی) و 
درنتیجه نزدیک شدن به جواب مطلوب است. در این پژوهش, ريشه 


عصلطم1)1 02 0110 -1 

(مصتعاونن هنق) ممتعافتاه اهنا -2 
۷ 01۷0)108 ۲۱۷ مهنص1 3-12 

طممص دععص]ز۳۲ -4 


نسبی مربعات خطا" (31۳5۳۴ل) به‌عنوان تابع برازش انتخاب شد. در 
گام چهارم» پس از انتخاب بهترین افراد از جمعیت حاضر برای نسل 
بعد» اصلاحاتی (تکثیر و جهش) روی جمعیت انتخب شده صورت 
می‌گیرد. در گام پنجم. فرزندان جدید تحت فرآیندهای پیشین در یک 
چرخه قرار می‌گیرند و برای تعداد معینی از نسل‌ها تکرار می‌شوند تا 
اه حل مناستبی بدا هنود (:() مر پیووهش خاضر از شرعفزر 
0 10018 060625010 برای به کارگیری مدل 0۳ استفاده شد. 
پارامترها و عملگرهای ژنتیکی مورد استفاده در محیط نرم‌افزار مذکور 
به همراه مقادیر آن‌ها در جدول ۴ ارائه شده است. 


توسعه مدل‌های هوشمند 


به‌منظور تعیین پارامترهای تأثیرگذار بر بار معلق " (م51) و تعریف 
الگوهاء از آزمون همبستگی پیرسون در محیط 5۳55 استفاده شد. 
مقادیر ضرایب همبستگی متغیرها در جدول ۴ ارائه شده است. نتایج 
حاکی از آن بود که کلیه متفیرها دارای تأثیر معنی‌دار در سطح یک 
۹ 
متغیرهای مساحت حوضه (۸) و زمان تمرکز (10) و هم‌چنین 
متغیرهای شیب )٩(‏ و ارتفاع متوسط (ظ) دارای همبستگی بالا 
(بزرگ‌تر از ۰/۷) بودند. بنابراین در تعریف الگوها از متفیرهای زمان 
تمرکز و ارتفاع صرف‌نظر گردید. نتایج آزمون همبستگی نشان داد که 
دبی جریان (0) بیشترین همبستگی را با بار ملق (لا) داشت. 
بنابراین به عنوان تنها متغیر ورودی در اولین الگو برای هر سه مدل 
هوشمند به کار گرفته شد. سایر متفیرها نیز بر اساس ضریب 
همبستگی اولویت‌بندی شده و در الگوهای بعدی به‌کار گرفته شدند. 
بنابراین بر اساس نتایج آزمون همبستگی, پنج الگوی مختلف برای 
تخمین بار معلق در گروه‌های اول و دوم مطابق جدول ۵ تعربف شد. 
شاخص‌های ارزیابی مدل 
به‌منظور آزمون عملکرد مدل‌های مورد استفاده در این مطالعه از 
شاخص آماری شامل ضریب تعیین (187), ريشه میانگین مربعات خطا 
(۷/5۲) و شاخص کارایی نش-ساتکلیف (5) مطابق با روابط زیر 
استفاده گردید (۴). 


۱۳( زور6 (و5-وک) ی 12 


2( ک) ۷ ۱۳ 


1 
۱۳ 2(بر5 - ن5)ا ۴ ]115 


02کصیا] زر را 
(۱۴) 1 > ۸۷ > هه ات | 1 - ۱5 


رامع مه معهتا50 100]16121۷6 -5 
س1) 1۵20 590۵2060 -6 


کاربرد سه مدل هوشمند در بر آورد بار معلق حوضه‌های آبخیز ... ۸۳۵ 


جدول ۳- مقادیر پارامترهای مورد استفاده در روش برنامه‌ریزی بیان زن 
۵00 6۱۳ م6 رز 0عون کامامهصن 200 وم)مصه و م6 ۵۲ ومباله۷ 11۳6 -3 م1201 


عملگرهای زنتیکی تنظیمات کلی 
63 ناه م6۵ 
رم 3 ۳۹ 5 ۱ ۰ ها 
0044 سر رت 30 تعداد کروموزوم 
مه ۵)20100 ]۷ 0۶ ۵۲طاصصناآ 
01 نرخ وارون‌سازی ر اندازه راس 
6 ۳۷۵5109[ 6 ۲1620 
01 نرخ ترانهش درج متوالی 3 تعداد ژن در هر کروموزوم 
۱۳۱ عم ۶ه مامتا( 
ممنانوممو112 ۴]5 یر ار تاره زره ۱۱۱۱۰۱۱۵۵ 
0۱3 نرخ ترکیب تک نقطه‌ای مجموع (+) تابع پیوند 
صمتاحصتططممع۲ اصلمموعو0 (ب) ممتان0ظ۸ صمنامص عصت‌اص[ 


جدول 2- مقادیر ضرایب همبستگی متغیرهای مورد استفاده در تحقیق حاضر 
۲ :۱ صز 0ععیه مماماهتنه۲۷ عط) ۵۶ اصمت‌لمی صمتاداهری ]0 معباه۲ 1۳06 -4 م1201 


بل <) 5 .1 3 5 ۸ 0 
1 0 
1 ۸ 
1 0283 ,5 
1 تک *ووفم ‏ 0214 ظ 
 06۵6 1‏ 9ق08 . تکف0 .و۲ 
1 0860 ت6کک ت06 02 5 
1 وا ا ‏ ا 0 ۱ ۵5(" 
1 0066 میم وا م۵0 ۵068 00 048 له 


2-021160(۰) 16۷۵1 0.01 عطا 26 عصهم‌تلتصونه مد ممتاماه‌ممن ۲ 


جدول ۵- الگوهای تعریف شده در تخمین بار معلق برای هر سه مدل هوشمند 
۵06۱6 )حهع1۳۱)61۱ معط معط تم 0و۱ 060صممرفیای عصت‌مهتای صا عمتتعصیء 0عصتلع0 116 -5 12016 


متغیر خروجی متغیرهای ورودی شماره الکو 
هه منم مماطاهانهه تمرح 0 متججم 
51 0 1 
11 0.۸ 2 
رل 0,۸,5 3 
0,۸,٩, ۶ 1‏ 4 
1 ۳( 5 


در روابط فوق. آا تعداد داده‌هاء مقدار ٩0‏ مقدار رسوب معلق اندازه 
گیری شده و 5 مقدار شبیه‌سازی شده برای رسوب معلق توسط مدل 
در گام زمانی ز ام می‌باشد. همچنین 5 متوسط مقدار رسوب معلق 
اندازه‌گیری شده می‌باشد. 


واسنجی مدل‌های هو شمند 
در پژوهش حاضر به‌منظور دستیابی به بهترین ساختار مدل شبکه 
عصبی, توابع انتقال (شامل تانژانت سیگموتید لگاریتم سیگموتید" و 


1-۸20 


خطی ) و تعداد نرون‌های لایه پنهان مورد ارزیابی قرار گرفت. طبق 
نتایج» برای هر دو گروه» تابع انتقال 1۸810 از خطای کمتری در 
مدل‌سازی برخوردار بود. عملکرد مدل ۸ به‌ازای الگوهای مختلف 
در گروه‌های اول و دوم و با استفاده از تابع انتقال مذکور در جداول ۶ 
و ۷ خلاصه شده است. مشاهده می‌شود که در هر دو گروه. به‌ازای 
اولین الگو عملکرد مدل بسیار ضعیف است و این موضوع حاکی از آن 
است که صرفا با بهکارگیری دبی جریان نمی‌توان برآورد صحیحی از 
بار معلق حوضه داشت. 


60( -2 
مصنله۲۳۷ -3 
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جدول "- عملکرد مدل ۸ به‌ازای الگوهای ورودی مختلف در گروه اول 
میم )وا ما رز مممت نمی مرح میمتوه م6 ما ۸۵ 66 0۶ ععصمصهمزنه۳۵ -6 12016 


18۶ ۰. 1۳61 )0۵/027( 5 


0.20 916 0.22 
038 9.42 031 
045 924 036 
040 906 0.40 
033 39 0۹4 


تعداد نرون الگو 
عصمنمه ۵ کمماصصت ۲ متقمعگ 
8 1 
7 2 
5 3 
6 4 
6 5 


جدول ۷- عملکرد مدل ۸۸ به‌ازای الگوهای ورودی مختلف در گروه دوم 
مرمع 60۵0و عظ) رز مومت ممتطاصی مرح میمتیو ما6 0۲ ۸۵۵۲ 6 0۶ رهم۳۵ -7 12016 


18۶ ۰. 1۳61۲ )0۵/027( 5 


018 ارطرا۱۱۱ 032 
023 10۳0601931 0۳46 
0.30 1(۳۹91 033 
032 10۳1909 نز 
034 10۳۹6.20 039 


طبق نتایج» با اضافه شدن مساحت حوضه به ترکیب ورودی» 
عملکرد مدل به‌طور قابل توجهی بهبود یافته است. به‌عنوان مثال در 
بپبود یافته است. به‌طور مشابه» در گروه دوم مقدار *12 در الگوی ۲ 
نسبت به الگوی ۱ از ۰/۲۲ به ۰/۴۶ افزايش يافته است. اضافه شدن 
سایر خصوصیات فیزیوگرافی حوضه به ترکیب ورودی‌ها نیزه به‌طور 
محسوسی تتایج مدل ۸ را بهبود داده است و در هر دو گروه 
بهترین نتیجة مدل مربوط به حالتی است که در ورودی مدل از کلیه 
متغیرها استفاده شده است. بهترین عملکرد مدل ۸٩‏ در گروه اول» 
مربوط الگوی ۵ با بیشترین ضریب تعیین (۰/۵۲-*10)» کمترین مقدار 
ريشه میانگین مربعات خطا (۷]۹۴-۸/۵۹۱0۵/0۵7) و بیشترین 
۵ بهترین نتیحه را به همراه داشت و مقادیر 122 :391 و 5 به- 
ترتیب برابر با ۱ 9 ۰۳۴ به‌دست آمتله ت برای 
گروه‌های اول و دوم تعداد نرون‌های مورد نیاز برای لایه پنهان به- 
ازای الگوی بهینه, به‌ترتیب برابر با ۶ و ۴ نرون بود. 

در ادامه به واسنجی مدل ۸۳۲5 جهت پیش‌بینی بار معلق 
تقسیم‌بندی داده‌ها در مدل ۸۳5 از روش افراز شبکه استفاده شد. 
این روش مبتنی بر انتخاب نوع تابع عضویت (متلشی» ذوزنقه‌ای» 
گوسی» زنگوله‌ای و...) و تعداد توابع عضویت برای هر متغیر ورودی 
می‌باشد. در این پژوهش, توابع عضویت مثلثی, ذوزنقه‌ای و گوسی به 
همراه تعداد مختلف آن‌ها مورد ارزیابی قرار گرفت و طبق نتایج به 


تعداد نرون الگو 
عصمنعص ۵ کمصاصصت ۲ متقععگ 
9 1 
6 2 
7 3 
5 4 
4 5 


دست آمده در هر دو گروه بهترین تابع عضویت از نوع مثلشی بود. 
عملکرد مدل ۸۱۷۲715 به‌ازای الگوهای مختلف در گروه‌های اول و 
دوم و با استفاده از تابع عضویت مثلشی در جداول ۸ و ٩‏ آورده شسده 
است. نتایج کلی به‌دست آمده توسط مدل ۸۳15 مشابه مدل 
۸ بوده و مشاهده می‌شود که اضافه شسدن خصوصیات 
فیزیوگرافی زیرحوضه‌ها به ترکیب ورودی‌هاء آموزش مدل را بهتر 
کرده و بهبود نتایج را به دنبال داشته است. بهترین عملکرد مدل 
15 در گروه اول مربوط الگوی ۵ با بیشترین ضریب تعیین 
(4)1۳-۰/۵۷ کمترین مقدار ريشه میانگین مربعات خطا 
(۹۳-۸/۷۲)۵۵/0۵]) و بیشترین ضریب نش-ساتکلیف 
(۱5-۰/۳۲) می‌باشد. در گروه دوم نیز الگوی ۵ بهترین نتیجه را به 
همراه داشت و مقادیر 127 ]۷۹ و 5[ به‌ترتیب برابر با ۰/۶۵ 
۷ و ۰/۴۰ به‌دست آمد. لازم به ذکر است که تعداد توابع 
عضویت برای متغیرهای 0 ۸۸ ٩:۰‏ و آطل) به‌ازای الگوی ۵ در 
گروه اول به‌ترتیب تعداد ۴ ۳ ۲ ۲ و ۲ و در الگوی ۵ گروه دوم برای 
متغیرهای مذکور به‌ترتیب تعداد ۴ ۲ ۰۲ ۴ و ۲ در نظر گرفته شده 
بود. 

همان گونه که قبلاً ذکر شد» یکی از گام‌های اساسی برای به 
کارگیری مدل 01:۳ انتخاب عملگرهای ریاضی برای تشکیل ساختار 
درختی است. بر اساس نتایج به‌دست آمده. در گروه اول به‌ازای 
الگوهای ۱ تا ۵ به‌کارگیری توابع مثلثاتی و جذر به همراه چهار عمل 
اصلی نتایج بهتری را به دنبال داشت, 
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جدول ۸- عملکرد مدل ۸۳۲5 به‌ازای الگوهای ورودی مختلف در گروه اول 
0۵ )کات معط صز مصمتامصتطاصصی )نامر متامتنه معط ۲۵ ۸۸۱۱۳۲5 6ظ) ۵ ععصمصصم] ۳۵ -8 12016 


1۳2 ۰ 1۳۷۲5۲ )۱0۵/027( 5 
0۳۵ 960 0026 
039 9.41 031 
048 913 038 
035 94 03 
057 813 05 


تعداد توابع عضویت الگو 
0 ص۱۲ متفمعگ 
4 1 
43 2 
3,2 4 3 
3,2,3 ,4 4 
2 ,2 ,3 ,4 5 


جدول -٩‏ عملکرد مدل ۸۱۳5 به‌ازای الگوهای ورودی مختلف در گروه دوم 
تمه 0صمععو عط) رز مصممصتماصی )مرا متامتنو عظ) رم ۸۸۱۱۳۲۵ عظ) 0۶ م‌صمص۳۵:]0 9 12016 


5 (۱0۵/027) 1۳۷۲5۳۲ .۰ 1۳2 
0.20 1091.0 03 
031 ۱۱00۳۵ 0351 
038 103117 060 
039 ااسل ۱۱0۳ 062 
0.40 10۳0107 065 


تعداد توابع عضویت الگو 
0 ۱ 
4 1 
42 2 
422 3 
4,224 4 
4,3 ,2 ,2 ,4 5 


جدول ۱۰- عملکرد مدل 217۳ به‌ازای الگوهای ورودی مختلف در گروه اول 
0۵ اکتا معط صرز مصمتاحصتماجصی )مرا مناماته معط ده بای معط ۶ن میصممصم]/۳۵ -10 12016 


1۳2 ۰ ۳۷۲۵۲ )0۵/02«7( 5 


968 031 
919 039 
14 0352 
1.93 033 
1.09 033 


الگو 
متحصعگ 
036 1 
045 2 
061 3 
064 4 
068 5 


جدول ۱۱- عملکرد مدل 0۳۳ به‌ازای الگوهای ورودی مختلف در گروه دوم 
0۵۵ 566000 66 صز مصمتمصنصاصصی )نامرصا منامتنه معط ها ای عمط ۶ن ممصهصمزن۳۵ -11 م1201 


1۷۲۲ )0۵8/02۲( 5 


1(396 028 
1۱۳1۰39 038 
1۳324128 041 
9382 042 
95.26 043 


در گروه دوم نیز ترکیب توابع لگاریتمی, نمایی» جذر و توان به 


همراه عهلگرهای اصلی با تنایخ خوبی همراه بود: عملکرد مدل طقی 
به‌ازای الگوهای مختلف در گروه‌های اول و دوم در جداول ۱۰ و ۱۱ 
نمایش داده شده است. همانند مدل‌های ۸۸۲۲ و ۸۳۲5 در مدل 
۴نآن نیز بهترین نتیجه مربوط به الگوی ۵ بود. به‌طوری که در گروه 
اول الگوی ۵ با بیشترین ضریب تعیین (418-۰/۶۸ کمترین مقدار 
ريشه میانگین مربعات خطا («00/0] ۳۹۳-۷/۶۹) و بیشترین 
ضریب نش -ساتکلیف (8-۰/۵۵) بهترین عملکرد را داشته است. در 
گروه دوم نیز به‌طور مشابه» الگوی ۵ به‌ازای مقادیر 57 1۷5 و 
5 به‌ترتیب برابر با ۰/۷۲ ٩۷۵/۲۶‏ و ۰/۴۳ بهترین نتیحه را نسبت 


1٩2‏ الکو 
۱ 
004 1 
056 2 
05 3 
09 4 
072 5 


به سایر الگوها به همراه داشت. 


مقایسه نتایج مدل‌ها 

با مقایسه شاخص‌های آماری عملکرد مدل‌ها به‌ازای بهترین 
سناریو (الگوی ۵) می‌توان دید که در گروه اول. عملکرد مدل ۳ با 
درا بودن بیشترین مقدار ضریب تعیین (۰/۶۸-*8), کمترین مقدار 
ريشه میانگین مربعات خطا (0۵/01) ۷]۹۳-۷/۶۹) و بیشترین 
ضریب نش-ساتکلیف (8-۰/۵۵) در مقایسه با سایر مدل‌ها پهتر 
بوده است. در خصوص گروه دوم نیز برتری مدل 013۲ با دارا بودن 
مقادیر 127 ]1۷ و 5[ به‌ترتیب برابر با ۰/۷۲ ۹۷۵/۲۶ و ۰/۴۳ 
مشهود است. 
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شکل ۵- نمودارهای پراکندگی مقادیر بار معلق مشاهداتی و پیش‌بینی شده توسط مدل‌های هوشمند به‌ازای بهترین الگو 
(الف)» (ب) و (ج): گروه اول (د)» (۰) و (وا: گروه دوم 
۵0 ۱0۵۵۵ 60۵ ۲0۲ ولم0مصظ )صمعناام)۱ عصمت ۱۵20 0ع0ه۵مرکیای ۵ عمباه۲ ۵60۱0660 0ج 0ع۲۷عموهان عظ) ۵ عاماص م5 -5 ۲1۵0۲۵ 
نامع 560080 1۳6 :) ورع) و(0) منامام افتت م1۳ <(0) و(رظ) و(ه) 


در شکل ۵ نمودارهای پراکندگی مقادیر بار معطق مشاهداتی و 
پیش‌بینی شده توسط مدل‌های هوشمند به‌ازای بهترین الگو در گروه 
در گروه اول. به‌ازای مقادیر کمتر از ۲۰ تن در روز, نقاط حول محور 
یک به یک متمرکز شده‌اند اما به‌ازای مقادیر بیشتر از ۲۰ تن در روز 
(تقریباً ۴ درصد داده‌های دوره تست) نقاط از محور یک به یک 
خاضله گر و حالی کیو امرخ ضل‌ها انار وود ]رس دا 
مربوط به مدل ۴ می‌توان دید که نقاط به نیمساز ۳۵ نزدیک‌تر 
هستند که نشان از تطابق بهتر مقادیر مشاهداتی و پیش‌بینی شده 
توسط مدل مذکور را دارد. در گروه دوم نیز به‌ازای مقادیر کمتر از 
۰ تن در روز تمرکز نقاط حول محور یک به یک دیده می‌شود و 
به‌ازای مقادیر بیشتر از آن (تقریبا ۱۲ درصد داده‌های دوره تست) 
پراکندگی ضعیفی مشاهده می‌شود و مدل‌ها حالت کم‌برآورد از خود 
نشان می‌دهند. با این وجود. پراکندگی و تمرکز نقاط حول محور یک 
به یک در مدل 0۳۳ نسبت به سایر مدل‌ها تا حدودی بهتر است. 
لازم به توضیح است که در هر دو گروه پایین بودن شاخص کارایی 
نش- ساتکلیف به علت پراکندگی نامناسب نقاط حول محور یک به 
یک می‌باشد. 


۱ 


«1 
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اعتبارسنجی مدل‌های هوشمند 

مدل‌های هوشمند در برآورد بار معلق بر اساس محدوده ضرایب تعیین 
() و نش- ساتکلیف (15") ارائه شده است. شکل‌های ۶ و ۷ میزان 
کارایی مدل‌های مورد استفاده در تحقیق حاضر را در تخمین بار معلق 
نمایش می‌دهد. با توجه یه شکل‌های مذکور, در گروه اول مدل‌های 
۸ و ۸۳۲5 با دارا بودن کارایی مشابه» از درجه دقت قابل قبول 
و رضایت‌بخش برخوردار هستند و به عبارت دیگر نتایج حاصل از آن 
بیشتری در تخمین بار معلق برخوردار است. در گروه دوم عملکرد 
مدل‌های ((۸ و ۸۳15 غیرقابل قبول است. در واقع با اینکه 
این دو مدل دارای ضرایب 1*7 خوبی هستند امابه علت دارا بودن 
مقادیربسیار پایین ضریب نش- ساتکلیف» نمی‌توان نتایج حاصل از 
آن‌ها را معتبر دانست. در مقابل» مدل طت1) از درجه اعتبار کافی برای 
تخمین بار معلتی گروه دوم برخوردر است. 
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شکل 7- میزان کارایی مدل‌ها در تخمین بار معلق برای گروه اول بر اساس ضرایب تعیین (122) و نش -ساتکلیف (5) 
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شکل ۷- میزان کارایی مدل‌ها در تخمین بار معلق برای گروه دوم بر اساس ضرایب تعیین (182) و نش-ساتکلیف (۱۷5) 
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برآورد میانگین بار معلق سالانه حوضه دره‌رود 

همان‌طور که قبلاً اثاره شد. یکی از مزایای مدل 0۳۳ ارائه 
معادله ریاضی بر اساس متغیرهای ورودی و خروجی می‌باشد. روابط 
۵ و ۱۶ معادلات به‌دست آمده توسط مدل 010۲ برای تخمین بار 
معلق را برای گروه‌های اول و دوم نشان می‌دهد. در معلالات مذکور, 
واحد متفیرهای دبی جریان (400 مساحت (۸ شیب (5) و بار معلق 
(:0)) به‌ترتیب متر مکعب در انیه. کیلومتر مربع. درصد و تن در روز و 
ضریب شکل (/5) بدون بعد می‌باشد. 


یی لک 
۳ تک ۶ ِ_ سب 
٩+ 510)-3.337+6( ِ (4 ِ ِ‏ 0 < 4 
(م5 < 2 + ۹/5 
- (6231 - )دب ها اند < ,ه 


۱)6( + ۱ ِ ِ | ۳" 


یکی از کاربره‌های عادلات :۱۵و ۱۶ که می‌توان آن‌هارا سل 
های منطقه‌ای رسوب حوضه دره‌رود نامید. برآورد میانگین بار معلق 
روزانه دبی جریان طی دوره قارع مشترک ۱۵ ساله (۱۳۸۰-۹۴) و 
خصوصیات فیزیوگرافی زیرحوضه‌هاء مقادیر بار معلق به‌صورت روزانه 
برای تمامی زیرحوضه‌ها محاسبه شد. سپس با تجمیع مقادیر بار معلق 
رو درا هر سییر سای تالف وتا با کین یر از 
مقا ال رسکیم پارسای ماه رما مسبت انم 
عاونا تیم هس کی بارش ماع ات 
زیرحوضه» دبی ویژه رسوب که نشان‌دهنده میزان فرسایش‌پذیری 
دره‌رود در جدول ۱۲ آورده شده است. با توجه به نتایج به‌دست آمده» 
سالانه حدود ۶/۳۲ میلیون تن رسوبات به‌صورت بار معلق در حوضه 
دره‌رود تولید و توسط رودخانه‌ها حمل می‌شوند. هم‌چنین میزان دبی 
ویژه رسوب در کل حوضه, ۱۰۴۷ تن در سال در کیلومتر مربع برآورد 
شده است. در تحقیق انجام گرفته توسط صمدزاده و همکاران (۲۱) 
مقدار رسوب برآورد شده برای حوضه دره‌رود. ۸/۷ میلیون تن در سال 
و میزان دبی ویژه رسوب ۶۲۴ تن در سال در کیلومتر مربع به‌دست 
آمد. تفاوت دوره‌های آماری مورد استفاده در تحقیق صمدزاده و 
همکاران (۲۱) (۱۳۵۶-۸۶) و پژوهش حاضر (۱۳۸۰-۹۴) می‌تواند 
تک از دلایان اخلاف ی تاه وس روش هرآ 
در تحقیق صمدزاده و همکاران (۲۱) از آمار بار معلق کلیه ایستگاه- 
های هیدرومتری حوضه دره‌رود استفاده شده است» در حالی‌که در 
پژوهش حاضر ماو تعدادی از ایستگاه‌ها (از جمله پل الماس 9 


گیلانده) برای مدل‌سازی به‌کار گرفته نشد. در سال‌های اخیر در 
فاصله‌ای نزدیک در بالادست ایستگاه‌های هیدرومتری مذکور بند 
انحرافی احداث شده است و این موضوع موجب عدم کاربرد داده‌های 
چنین ایستگاه‌هایی گردید. طبق نتایج جدول ۱۳ بیشترین حمل 
رسوب در سطح حوضه دره‌رود در ایستگاه بران واقع در رودخانه دره 
رود رخ می‌دهد و ایستگاه آتشگاه واقع در رودخانه نوران‌چای کمترین 
فرسایش ویژه را به خود اختصاص داده است. نکته قابل توجه این 
است که در سرشاخه‌های حوضه دره‌رود که شیب حوضه زیاد است. 
مقدار دبی ویژه رسوب کمتر بوده و هر چقدر به خروجی حوضه 
دیجم شویم و هرفن ۱ تکام »روم که ها شیب خومه ماهتن 
پیدا کرده و دبی ویژه رسوب بیشتر می‌شود. این موضوع قبلا در 
آزمون همبستگی نیز دیده شده بود و ضریب همبستگی بار معلق با 


شیب حوضه مقدار منفی به‌دست آمده بود. 


نتیجه گیری 

در این مطالعه به بررسی عملکرد مدل‌های ۸۵۳۲ ۸۱۳۲5 و 
۳ در پیش‌بینی میزان رسوبات معلق حوضه آبخیز دره‌رود وآقع در 
استان اردبیل پرداخته شد. طبق نتایج» ضعف مدل‌ها در تخمین بار 
سل ضرف با اسفاده از مقادیر هی هیام (الکوی ۱ تمابای بود و نا 
اضافه شدن خصوصیات فیزیوگرافی به ترکیب ورودی‌ها نتایج مدل‌ها 
بهبود یافت. که از این لحاظ با نتایج به‌دست آمده توسط ابوالفتحی و 
همکاران (۱) مطابقت دارد. در مقابل بهترین عملکرد مدل‌ها در 
تضمین بان ملق متفلی به الکو ۵,عامل ت کیب ورودی 0 5 
۰ و () بود. در گروه اول مقادیر شاخص‌های ۳2 ۳5 و 113161 
به‌ازای بهترین سناریو (الگوی ۵) برای مدل ۸ به‌ترتیب برابر 
و ۸۵ قن فر روز برای تال 018 رنه گرتشب پرایی 
۲ دقن موه ویزای معل ۵0۳ نتهم کیب پراب: 
۸ ۰/۵۵ و ۷/۶۹ تن در روز به‌دست آمد. هم‌چنین در گروه دوم 
مقلدیر شاخص‌های مذکور به‌ازای بهترین سناریو (الگوی ۵) برای 
مدل ۸ به‌ترتیب برابر ۰/۵۹ ۰/۲۴ و ۱۰۴۶/۲۰ تن در روز برای 
مدل ۸۳۲5 به‌ترتیب برابر ۰۰/۶۵ ۰/۴۰ و ۱۰۰۱/۰۷ تن در روز و 
برای مدل ۳۳) به‌ترتیب برابر ۰۰/۷۲ ۰/۴۲ و ٩۷۵/۲۶‏ تن در روز به 
دست آمد. بر اساس نتایج اعتبارسنجی و با در نظر گرفتن شاخص 
های ضریب تعیین و نش- ساتکلیف» کارایی مدل‌های ۸۷ 
۸5 و 0۳۳ در شبیه‌سازی بار معلق ایستگاه‌های گروه اول قابل 
قبول بود. اما در گروه دوم فقط نتایج مدل 0۳۳ معتبر بود و عملکرد 
مدل‌های ۸۲۸۲ و ۸۱۳15 غیر قابل قبول بود. 
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جدول ۱۳- مقادیر میانگین بار معلتق سالانه و دبی ویژه برای کل حوضه دره‌رود 
۵ ۱2۲۰۵۱۱۰-0 46 ۲۵۲ عمط من عل6هری 0صه 0م۱ 060جهمرعیاه آمیتصصه عمط ۵۶ و۷2۱6 -12 12016 


دبی ویژه رسوب میانگین بار معلق سالانه 
ملع 0 واه 0 ۸۲۵۲۵۵ 
ای /۱0۵) ععتمط‌عنل (۵۵/۲۵۵۲) 1020 
526.12 113.7 
5۹04 2000 
508 13920 
4 10 09.32 
32-1 212 
1۳05 22.02 
1285 199137 
926 12967 
9.72 1-933 
12359 109930 
6408 6.08 
1632 2-3039 
11003 62۹223905" 
193 113-93304 
19-14 1۱۳۱9380 
19373 2-2921 
٩62۰2323210۴ 10۳9۹01‏ 


در تحقیق حاضر مدل 0۳۳۴ بهترین عملکرد را در تخمین بار 
معلق حوضه به همراه داشت که در تحقیقات مشابه انجام شده توسط 
کیسی نو شیرین (۱۴) و یی پور و همکاران (۳۳) نیز این مدل 
بهترین نتیجه را مقایسه با مدل‌های ۸ و ۸۸۳15 نان داد. 
ضعف عمده مدل‌های هوشمند در برآورد مقادیر پیک بود و حالت کم 
برآورد مدل‌ها در تخمین مقادیر بیشینه بار ملق نمایان بود. در 
تحقیق مشابه انجام گرفته توسط کیسی و شیری (۱۶) نیز ضعف 
مدل‌های ۸ ۸۳۲5 و 0۳5۳ در پیش‌بینی مقاد 


معلق دیده شد. بر اساس الگوی بهینه با استفاده از مدل ۳۲۳ برای 


پیب مساحت نام ایستگاه نام رودخانه 
(60) 06م۱٩‏ (فصمل وم مصهه صمتاهاه عصوه 135۲7۵ 
1111 7373 ننه کران نرگس‌چای 
۱۳۱۵۱۹ ۱۱2865-2۷7 
7212 24975 تاه نوران چای 
طحعطءه‌اض / ۱۱ 
237 28 ان قره‌سو 
120170 0212-0 
5314 28112 عنبران عنبران 
و( و( 
۱ قزچی چا 
298 82319 بریل سقزچی‌چای 
11 «دط-نط520062 
2.64۱ 13122 هیر هیرچای 
۳11 «حل۲۲ 
5.71 13158 نوران نوران‌چای 
۱۹۹ -۱۲۱۲۵۵( 
ز خیاوجا 
دوم 140029 پل سلطانی یاوچای 
6 501201 ۷- ۱6۳۱۱3۷ 
۹833 16333 فیر نیرچای 
۲ «حل۲۱۳( 
5113 6 یامچی بالخلی چای 
تطم۲2 )بو تاهظ 
1031 804402 ۱ کوزه تبراقی قوری‌چای 
نهه0۵۳] 100761 0۵00-02 
4و ۰ 4066829 سامیان و 
صعتصع -0۵۲۵۳ 
تو ...هجو جر آریاب‌کندی قره‌سو 
رصق واوصا-ت۸ر -طفه0۵ 
0473 7681 دوست پیگلو قره‌سو 
100 10056 -طفنه0۵ 
9 ۰ ۰ 9600024 مشیران درهرود 
صح‌تنطوه]۱۷ 1270-00 
85 ۰ 12206744 بران دره‌رود 
۱ 1270-00 
مجموع 
او 


هریک از گروه‌های اول و دوم مدل‌های منطقه‌ای رسوب (معادلات 
رنه خر خی انم تفه فان سای سر گرد 
رود در رودخانه‌های تکامل يافتة آن یعنی دره‌رود و قره‌سو تولید و 
منتقل می‌شود. نظر به اینکه اندازه‌گیری بار معلق در ایستگاه‌های 
هیدرومتری» به‌ویژه در موآقع سیلابی با مسائل و مخاطراتی همراه 
بوده و به‌صورت منظم انجام نمی‌گیرد لذا از معادلات منطقه‌ای 
رسوب می‌توان برای تخمین رسوب معلق حوضه دره‌رود در روزهای 
فاقد آمار بهره برد. جهت تعمیم نتایج به‌دست آمده در این مطالعه. 
باید حالت کم‌برآورد مدل‌ها مورد توجه قرار گیرد و به‌منظور بهبود 


۲ شریه آب و خاک. جلد ۳۴. شماره ۴. مهر - آبان ۱۳۹۹ 


نتایج حاصل از تحقیق حاضر پیشنهاد می‌شود کارایی مدل‌های رودخانه‌ها تا حد زیادی متأثر از ذوب برف می‌باشد» انتظار می‌رود در 
این موضوع بررسی شود که پیش‌پردازش (تجزیه) داده‌ها توسط تبدیل این مسأله بر روی سری زمانی داده‌ها تأثیرگذار بوده و ممکن است بر 
موجک!" تا چه حد می‌تواند به بهتر شدن نتایج مدل‌ها کمک کند. از عملکرد مدل‌های ارائه شده نقش داشته باشد. لذا می‌توان چنین 
سوی دیگر با توجه به تأثیر تغییر اقلیم در کشور که با کاهش بارش و مطالعاتی را در دوره‌های معین چندین ساله. جهت به‌روز رسانی مدل 
افزايش دما همراه بوده و در حوضه آبخیز مورد مطالعه» دبی جریان با شرایط منطقه تکرار نمود. 
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0ص ,(۱۵۵/027) 1000 صقطا وفع ۵۶ ومببله۷ ما 10۶ صععو فه ممتوهمون اصماامعه مطا ,نامع فصممعو فظ) 0 0۵96۲۷۵0 
طان مط ۵۶ ممتوومعن مط روونامع طاهه ص ۳۱۵۱۷۵۷ .فصمتاتصمی صمتاهمدتاوم‌علهه 20 واممصظ فطا رلهتعصمع 1 
5مج ۸۲۳18 0صج ۸۵۲ ,کل 0ص2 2 ۶ه ومبلج* عمط من 2560 .عام0مصه رمطان عط صقطا تعتاه غقط«عصصمو و72 مج 
مره مه ماصادتاه: عمج مه اع0مصص نزن قط1 .منامتع )وت مطا ب10 مدمه بامامه‌گفتامه صه ماطامامعععه مطا 20 
۵۵ ۸۵۱۳۲۱۶ 40صد ۸۵۱۵۲ ۵۶ ب«مجمتمتلله معط ,صقط تعطام فطا مر .منامعع اف فطع ۵۶ 1۵20 0ع0صومویاه عط) عصتافطتاده ظ1 
اصص معط مطا عصلوته واعهمصه عصععع ۵۶ فاناوع: مطا گه مموتنمموصصمی .متامتع. فصصععو فص ع10 )20000 2۵۶ 
0۷ 10۷/65 عطا ,(0.68 - ۲2) امعت‌الهمن مملافصتصمامل )ومطونط مغ ۱1 ۵06و نان مطا عقطا 0عهع4ص1 ممتاقمتماصمع 
ماه ط«مطه فقط وماو صمئلهل ۷۵1 مطا صا 0.55 مه لقناهع ول م1 (رد/۵8 7.69 ح ۴9) مه عمفتان0و 32 
0 136 ۲۵۲ ,رتم۱ نصا .منامع افیا فطع عم 1۵20 60)صوم‌ویاه فطل عصتامصتاوه صا واع0مصه ععطم فط صقطا معصفممهم 
عم عتفناود مصفعط اموز ]109/6 عطا ب(0.72 ع 2) صعتم‌تااعمی صمتاعصتصعاعل )ععطونط عمط طفنه اعمصه نام مط روتاونع 
تعطاه صقطا عمصفصتهمگهم ععااهها مبومطو فقط وماه ممتاه۷211 مط ما 0.43 م6 لقتاوع ول 16 .(1۵۵/02 975.26 2 /۴(۱9]۲) 
0۰ 60 جووویاه عط عصتلهصصتاوی صا واع0‌مظظ 

معط صد 0مامعتاو1۳۷6 ۳/۵۵ فاعم اصمعتااماصد احمعهلتل زه تمممته‌تاله مطا رای تصمععتم فط ما خصمتعتمومن 
۶ (۷6۵۲۶ 15 وصتتا ریع.1) 0۵۲100 060ه۵ه صرح رلتهعه: ونط) ص .0مطومه ۲0ا0عجمنونا ۵۶ ممتاقصصتاهء 1020 060ظ۵م‌فناه 


طممطممطه۷۱ که لویع نا رمصتمممتمص ماج ۷۷ ۵۶ اصمصصندمعن1 روع80و۲۳۵16 اصفادولفه 20 ,۷۲,۹6 0ماهنهرن -3 20 2 ,1 
۰ ,11۳27 ,۸۳02011 ,۸02011 

(15. 26 .مصصن ۵ 1100 _ مه تاتممصط عمطمه مصت‌طمصدمتزمن -۶) 

1 طومطممطم۱۷۱ ۵۶ اهتنا رومعتهفمک آهتبمداط ی متبانه‌ه۸ ۶ معقاامن صقطع ۱۷۵ ,۲۲۵16990۶ مع001وو۸ -4 

1۵]: 1022067/5۷. 1 


کاربرد سه مدل هوشمند در بر آورد بار معلق حوضه‌های آبخیز ... ۸۴۵ 


۸۵ 3 0 1۵62060 فصمتاهاد متتاعطصمل بط عمط 2 )صمصتلعه فصه مععمطم‌ونل «م۶۱ مصتل‌بدامصد رقئعل 2۵۵8۲6۵0 
۵ 0 تاقصاطاصصم . انامرصا. امن وب رجم1 . 0ع0صمم‌واه. فطا متقاصته مه بعهتم ‏ 1 وا موه رلعطفماه۳ 
ع عصتمتطاصم ها 7۵0 متمرح قه۷ ممتاقصصتاوه ‏ 1۵20 0ع4هصومفناه ۵ پ«مجتتاهعه فطع رقام0مصظ له ۲۵۲ ,00۲60قص۵هم 
لمصمنعع1 رباعم نان معط ۵۶ بمدتنهمه تمطونط مطا مه بان ,عمییلهه موعقطه‌فیل مه دما فتهامم مه ملطامرهتومتوردام 
0 1020 460مممرویاه اقتتصصه ,موه .7 املههو۵و رومتامتع ۹۵6080 4صه ف ما 10 7۵0منطمج ۷۵۲۵ واععمصظ عصممصتعو 
0 عطا 0۶ ]005 رفالناوعع عط ما عصتلت0ععض .فصتفوه- ماه اه 10۲ 0متملنملوه مه معتهطمفنل ملموصرو اصممصت8۵0 
02121-50(۰) 0ج تام-طمت1۲2 ریعی) ۲۱۷۵۲9 010 168 صا 0ما۲ممفصها هه 0100660 15 ۷۸2۵۵۲8۴۵۵0 0بامک-طم۱2۲ عمط 0۶ 1020 
فصمصصتلمو 060مممرویاو ۵۶ قمما ممتااتص 6.33 رلمطوهماه منامع-ط۵ت2نا مطع طد رط‌تهعوع۲ ونطا ۵۶ فالناوع عطع مه ۳860 

2001-5۰ مطلتتال 0۲)60مفصهتا ۷۸۵۲۵ 


٩60116۶ 6۵‏ 9۵۵01610 رواعمظ )186۵11126 رلتامک-دامتهنا ر1۵2 0ع4هممفناو: :1۵۲۱۷۵۲۵6 


